Ecology and Evolution

The health and wellbeing of living organisms and the environment are irrevocably interconnected. Animals constantly shape the environments in which they live, and environmental factors, in turn, shape the minds and behaviours of animals. The emerging discipline of environmental neuroscience studies this reciprocal interaction and offers a unique opportunity to integrate sustainability into undergraduate training. Here, we present and evaluate Forging Sustainable Research, an interdisciplinary undergraduate course that integrates teaching of neuroscience, biology, psychology, environmental science, and sustainable development. Students explore how research in these fields can be interlinked 24 and re-oriented towards addressing emerging environmental challenges. In parallel, students learn how to measure and mitigate the environmental impact of research activities. Student feedback suggests that the course is effective at raising awareness and understanding of sustainability in research, while also enhancing student perception of their own ability to promote positive change. As climate and ecosystem change impact all sectors, sustainability awareness in academic and professional spaces is increasingly sought after. Students in the course therefore receive training in transferable skills relevant for diverse career paths. In the long term, this course may serve as model for direct integration of sustainability training into undergraduate curricula in the natural sciences.

Fluctuating asymmetry (FA) is widely used as a proxy for developmental instability, although reduced rather than increased FA, has occasionally been reported under particularly challenging conditions, possibly because environmental filtering removes less developmentally stable individuals. Here, we investigated temporal variation in wing shape asymmetry in males of Drosophila sturtevanti collected in 2022, 2024, and 2026 at the same sampling site, an urban Atlantic Forest fragment. Sampling in 2022 and 2026 were completed during single field expeditions in early April, whereas in 2024 repeated collections were required until mid-May because of markedly low capture success. We therefore examined whether this difference was associated with unusual environmental conditions and whether FA was higher in 2024, as expected under greater developmental perturbation. Daily temperature and relative humidity were examined from early March onward, encompassing approximately one month before adult emergence. The 2024 sampling period was warmer and drier than the corresponding periods in 2022 and 2026 and coincided with the final phase of the strong 2023-2024 El Ni&ntildeo event. Contrary to our prediction, FA was significantly lower in 2024. The magnitude of directional asymmetry remained similar among years, although the spatial pattern of the DA vector differed between 2024 and 2026. The coincidence of reduced capture success and reduced FA is consistent with the hypothesis that adverse conditions may impose demographic filtering, such that individuals expressing greater developmental instability are less likely to contribute to the sampled adult population.

Rising temperatures linked to global climate change can disrupt the chemical communication that animals rely on for reproduction. Arthropods are particularly vulnerable due to their limited ability for thermoregulation and their reliance on cuticular hydrocarbons (CHCs), which have a dual role in social communication and heat protection. Heat-induced changes in cuticular composition may therefore lead to mate recognition errors. We tested this in the burying beetle Nicrophorus vespilloides, by exposing males to either benign temperatures (20{degrees}C) or a simulated 3-day heatwave (26{degrees}C). We recorded the frequency of same-sex sexual behaviour under these contrasting thermal conditions and examined associated changes in the males' cuticular profiles. Our results show a higher rate of same-sex mounting under heat stress, with the CHCs of heat-stressed males shifting toward more saturated, longer-chained compounds and more female-like compounds. We also found a positive association between chain length and the frequency of same-sex mounting, suggesting a trade-off between the social communication and waterproofing functions of CHCs. Overall, our findings provide evidence for a causal effect of heat on same-sex sexual behaviour, which may be partly explained by changes in cuticular lipids. Such disruption of chemical communication could have major consequences for animal reproduction in a warming world.

When mating among genetically related individuals results in inbreeding depression, various behaviors and life-history traits evolve to favor outbreeding. Long-range sexual communication following natal dispersal is one such adaptation. In parasitoid wasps with single-locus complementary sex determination (sl-CSD), inbreeding increases the production of unviable or sterile diploid males. Mastrus ridens is such a species, with a genetic load underpinned by sl-CSD, combined with a high potential for sibmating because siblings emerge in close proximity from the same host and become sexually active shortly after emergence. With this context in mind, we investigated the mate-finding strategy of M. ridens using complementary behavioral, chemical, and electrophysiological approaches. Wind-tunnel bioassays and field experiments showed that males orient toward virgin females and, to a lesser extent, mated females. This attraction from a distance appeared to be the sole form of conspecific attraction. Chemical analyses using static and dynamic headspace sampling and solvent extracts, followed by GC-MS, revealed several potential pheromonal candidates, including methyl 6-methylsalicylate. This compound elicited reproducible antennal responses in males in GC-EAD recordings. Altogether, our results suggest that mate finding in M. ridens is guided by a volatile sex pheromone emitted by females, with methyl 6-methylsalicylate as one putative component. Such mate-finding mechanisms promote outbreeding and may explain the low proportion of unfit diploid males observed in natural populations. Surprisingly, chemical analyses suggested that males also emit this compound, whereas neither males nor females showed behavioral responses to males. This apparent discrepancy suggests that attraction may depend on a more complex pheromonal signal and opens further avenues for investigation.

Riverine and riparian food webs are connected via 'allochthonous' resource flows from neighboring ecosystems. These resources, such as terrestrial leaf litter, contribute to benthic macroinvertebrate diets, supporting higher trophic levels. While allochthony occurs in streams globally, its sensitivity to broad environmental gradients remains poorly understood. We explored the drivers of allochthony across six river basins within the contiguous United States, including climate, stream size, and flow permanence. We used Bayesian mixing models to estimate macroinvertebrate dietary proportions in four functional feeding groups. Air temperature best predicted allochthony, with consumers relying less on allochthonous resources with increasing temperature. Our models indicate that macroinvertebrate diets become primarily autochthonous at approximately 15{degrees}C mean annual air temperature, with leaf litter-consuming shredders showing a slightly higher threshold for this shift. We also found increased allochthony in smaller, non-perennial streams. These findings highlight the complexity of stream allochthony and responses to multi-scale environmental factors.

Social aggregation is widespread in vertebrate larvae, yet how genetic relatedness shapes group cohesion remains poorly understood. In particular, the association between genetic similarity and social aggregation may depend on the scale at which genetic relatedness is evaluated. Here, we quantitatively examined aggregation behaviour in tadpoles from three lineages of Japanese toads, Bufo japonicus, B. formosus, and B. gargarizans miyakonis . We experimentally combined tadpole cohorts differing in relatedness and lineage and quantified aggregation using spatial dispersion analyses. We then examined relationships between aggregation behaviour and genetic variation using mitochondrial DNA (mtDNA), major histocompatibility complex (MHC) class II variation, and microsatellite diversity. Overall degree of aggregation did not differ significantly between sibling and non-sibling cohorts, and Morisitas index provided no evidence that tadpoles formed separate clumps according to relatedness. However, the correlation between genetic variation and aggregation differed across scales of genetic relatedness: among non-sibling cohorts, greater mtDNA pairwise distance was associated with weaker aggregation, whereas among sibling cohorts, greater MHC amino-acid difference was associated with weaker aggregation. These results indicate that the degree of aggregation is related to genetic similarity across different scales of relatedness, rather than being explained by a simple sibling/non-sibling distinction. Our findings suggest that tadpole group cohesion is related to genetic similarity in a scale-dependent manner, highlighting the importance of considering the scale of genetic relatedness when investigating social discrimination in vertebrates, while the sensory and molecular mechanisms underlying these relationships remain to be determined.

The rapid global expansion of dams has significantly reshaped freshwater ecosystems, yet long-term trophic dynamics and food-web stability within reservoirs remain underexplored. This study examines the reorganization of fish communities induced by the trophic surge following dam construction and river impoundment. We used the predator-to-consumer mass ratio (PCMR) as a proxy for trophic pyramid structure to track changes in fish biomass distribution across three large boreal reservoirs in Quebec, Canada, over 20 years, using data from before and after impoundment. We found transient top-heaviness events (i.e., PCMR values > 1; predator accumulation), characterized by a hump-shaped trajectory of the PCMR values in two reservoirs (Opinaca and Robert Bourassa). The third reservoir (Caniapiscau) showed a gradual increase in top-heaviness but did not exhibit a clear transient pattern. The magnitude and duration of top-heaviness varied across reservoirs, potentially due to differences in the system's capacity to receive and accumulate nutrients, which are driven by water residence time and reservoir size. This research provides a conceptual framework for understanding the long-term trophic dynamics of reservoirs, which can help inform the development of more effective management strategies to preserve ecological health in these ecosystems.

Forecasting infectious disease dynamics is important for ecological management and emerging disease preparedness, yet many systems lack the long-term datasets required to develop reliable predictions. Spatial replication may provide an alternative by allowing information to be shared among populations, but the extent to which this improves forecasting remains unclear. Here, I used epidemic and temperature data from 20 replicated semi-natural Daphnia-parasite pond populations monitored across four seasons (80 epidemics) to test whether information from multiple populations improves forecasts and whether environmental similarity identifies informative populations. I forecasted disease prevalence, infected host density, and healthy host density using a benchmark model, autoregressive integrated moving average (ARIMA), and time-series regression models. Models were trained using focal-population data alone, mean dynamics across other populations, or temperature-weighted mean dynamics from environmentally similar populations. Epidemics showed strong seasonal structure, with prevalence and infected host density generally peaking during the warmest period of the season. However, no forecasting approach consistently outperformed others across all variables. Temperature improved forecasts of prevalence and infected host density but not healthy host density. Multi-population training improved forecasts in some cases, particularly for ARIMA models, whereas temperature-weighted averaging provided no additional benefit over simple averaging. Regression models consistently produced the most accurate forecasts of healthy host density. These results show that borrowing information among populations can improve forecasts under specific conditions but is not universally beneficial, highlighting the importance of matching forecasting approaches and data sources to the ecological processes being predicted.

Precocious reproduction, prior to the acquisition of the final adult morphology, occurs in a variety of animals and represents a distinct life-history strategy with important consequences for growth and survival. Larval reproduction is far rarer and usually restricted to groups that go through metamorphosis. In ctenophores, Carl Chun described dissogeny more than 130 years ago, a distinctive life-history pattern involving larval and adult reproductive periods separated by gonadal regression. Yet, the potential plasticity of larval reproduction in ctenophores remains unresolved. Here, we developed a method to follow and quantify reproductive output in Mnemiopsis leidyi larvae and tested environmental factors that may regulate precocious reproduction. We found that conspecific density strongly influences its onset: low-density, well-fed animals can bypass reproduction at the cydippid stage and develop rapidly into the lobed form, whereas higher density promotes precocious reproduction. Food availability and associated growth rate also contribute but cannot alone explain the observed differences. Once reproduction is established, however, high conspecific density suppresses reproductive output, an effect reproduced by exposure to seawater incubated with high M. leidyi densities, consistent with water-borne chemical communication. Salinity and light conditions further modulate spawning. Together, our results reveal substantial environmental plasticity in precocious reproduction in M. leidyi, reconcile long-standing contradictory observations of ctenophore development, and highlight how environmental cues can reshape reproductive life-history transitions.

Immune defenses may vary across environments due to differences in life history or abiotic and biotic characteristics, but few studies have compared the strength of immune responses across broad geographical gradients. Here we investigated how bacteria killing ability (BKA), a measure of constitutive innate immune function, varies among populations of tree swallows (Tachycineta bicolor) inhabiting different environments. We hypothesized that populations living in unpredictable environments, which also have shorter breeding seasons and thus fewer opportunities to renest in the event of reproductive failure, would invest in reproduction at the expense of immunity. To test this hypothesis, we measured BKA and monitored reproductive activity in four populations located across most of the breeding range of tree swallows, in Tennessee, New York, Wyoming and Alaska. Consistent with our predictions, BKA was lowest in Alaska and Wyoming and highest in Tennessee. In most populations BKA was also lower during incubation than during nestling provisioning, and lower in males than in females. These results suggest that the high reproductive investment seen in populations experiencing more unpredictable breeding conditions may come at a cost to immune function. While reduced immune investment could increase susceptibility to parasites, this potential cost is likely offset by lower parasite pressure in these environments, limiting the need for maintaining strong constitutive innate defenses like BKA.

Vertebrates cope with rapid temperature fluctuations by plastically adjusting their behaviour and physiology. This phenotypic plasticity is partly mediated by plastic changes in glucocorticoid levels; however, whether individual differences in endocrine plasticity are associated with fitness remains unknown. Here, we examined the relationship between individual variation in glucocorticoid plasticity in response to air temperature changes and reproductive success in a wild great tit population (Parus major). Over nine years, we collected repeated measurements of baseline and stress-induced corticosterone during nestling rearing. We found consistent reaction norm patterns in hormone response to temperature variation -- as shown in a previous study conducted over a shorter period. This provides robust evidence of individual differences in average hormone levels (intercepts) and the strength of the plastic response (slopes), indicating potential for evolution. However, neither intercepts nor slopes were associated with reproductive success (fledgling number or brood mass), suggesting that individual variation in endocrine plasticity was not under natural selection, at least within the environmental conditions experienced. This study provides a first empirical test of the evolutionary implications of endocrine plasticity in the wild and highlights the need to assess its fitness relevance under more extreme conditions, such as those expected under future climate scenarios.

Whole genome analysis of spittlebug populations, an Aphrophora species (Hemiptera: Aphrophoridae) known to spread Pierce's disease, from various vineyards in Napa and Sonoma Counties demonstrates that these insects are conspecific with Aphrophora permutata found in quasi-natural populations. The genomes analyzed fall in three autosomal clusters with two sympatric populations (based on admixture) that diverged around 40,000-90,000 years before present according to demographic models. We also recover three related but distinct mitochondrial clusters. There is no genetic differentiation between Napa and Sonoma County vineyard populations and very modest differentiation among vineyards and quasi-natural Berkeley/Oakland Hills population. These vineyard and quasi-natural populations have distinct life histories, including divergent oviposition hosts (understory herbs versus pines) and different patterns of adult seasonal migration (out-migration and return versus simple movement from the understory to the pine canopy). Analysis of the mitochondrial gene COI indicates that A. permutata, Aphrophora fulva and Aphrophora maculosa are intermingled in the same mitochondrial clade, along with the more distantly related species Aphrophora gelida. Further, all three autosomal genetic clusters show evidence of selection for insecticide resistance. We also find that endosymbionts Wolbachia and Rickettsia were both present in a substantial fraction of specimens analyzed, but neither appears to influence population structure. Overall, our findings elucidate an intriguing degree of life history divergence between geographically proximate, conspecific populations linked by gene flow and may have implications for the continued management of vineyard populations.

1. Soil invertebrate biomass is closely linked to ecosystem processes, yet its routine quantification remains constrained by the effort required to weigh large numbers of small individuals. Length-mass relationships provide an indirect alternative, but body length can be difficult to standardise and extract automatically from images of curved, coiled or contracted organisms. Projected body area may provide a more tractable image-derived predictor, but its performance across diverse soil-faunal taxa and at the community scale remains insufficiently evaluated. 2. Here, we developed a framework using projected body area and contour length to estimate fresh and dry mass from measurements of 1,792 soil invertebrates representing 15 higher-level taxa. We fitted taxon-specific models for six well-represented groups and general models across taxa, and evaluated predictive performance using leave-one-out cross-validation. We then assessed cross-dataset transferability by comparing our data with those of Lemoine et al. (2026) and pooled the two datasets to derive combined models. Finally, we tested whether individual predictions could be scaled to independently measured community biomass. 3. Taxon-specific models generally reduced prediction error relative to general models, with mean absolute percentage error (MAPE) values of 6.9-37.3% for fresh mass and 7.6-30.2% for dry mass. General models retained MAPEs of 27.2% and 22.4%, respectively. Comparisons of our dataset with that of Lemoine et al. (2026) showed no detectable differences in six of eight taxon-by-mass slopes. Although intercepts frequently differed, no overall dataset effect was detected after accounting for taxa. Pooling both datasets broadened model development, while taxon-specific models generally retained lower error than the combined general model. Finally, community-level predictions tracked measured biomass across all model sets (R2 = 0.87-0.99), but absolute error varied among models and faunal groups (MAPE = 9.6-52.5%). 4. Body area and contour length extracted from images therefore provide a practical basis for scaling soil invertebrate biomass estimation from individuals to communities. Integration with automated imaging and segmentation tools such as flatbug, BIODISCOVER and CollembolAI could extend soil invertebrate monitoring at larger scales beyond abundance and taxonomic composition to biomass, while providing scalable inputs for analyses of food webs and general ecosystem ecology.

Self-fertilization results in reduction of heterozygosity over time, accelerated by small population sizes, restricted ranges, and fragmented habitats. Yet some species maintain heterozygosity despite these factors. One is Corallorhiza bentleyi, a threatened, mycoheterotrophic, self-pollinating orchid restricted to Virginia and West Virginia, USA. Reduced representation sequencing was employed across five of the ~12 known localities of C. bentleyi, including representatives of sister species C. involuta from Mexico, to quantify patterns of population structure and genomic diversity. Surveys of C. bentleyi revealed a predominance of the cleistogamous (closed-flowered) growth form. Differentiation was observed between the two species but with weak population structure in C. bentleyi. Surprisingly, excess heterozygosity and negative inbreeding were recovered in both species, counter to expectations. Analyses provided evidence against bioinformatic artifacts (collapsed paralogs), small-sample bias (Robertson effect), and accumulation of mutations along chromosomes due to low recombination associated with selfing (Meselson effect) as explanations for excess heterozygosity. Demographic analyses favored a model of post-glacial divergence among C. bentleyi and C. involuta lacking bottlenecks. Finally, simulations under alternative selection scenarios favored models of overdominance vs. those of negative frequency dependent selection. It was concluded that overdominance, namely selection against selfed homozygotes, is the most likely explanation for elevated heterozygosity in C. bentleyi. Conservation concerns are discussed in light of these findings and the threats of extinction for one of the most vulnerable orchid species in North America.

Researcher disturbance can alter the behavior, survival, and breeding success of wildlife, with effects that are especially well-documented in waterbirds. Anecdotal reports suggest that disturbance can even cause cavity-nesting species to alter the structure of their nests. We used a long-term study of a burrowing seabird, the Leach's Storm-Petrel (Hydrobates leucorhous), to ask whether storm-petrels renovate existing burrows and whether such renovation is related to a history of researcher disturbance. Our study was conducted at the Bowdoin Scientific Station on Kent Island, New Brunswick, Canada, where in 2024 we measured the burrow architecture of 363 burrows, including length and structure for occupied and unoccupied burrows both inside and outside a long-term study site with consistent research activity. Of these, 150 burrows had been measured in 2006, allowing for direct longitudinal comparison of burrow architecture. Nearly all burrows changed in length over the 18-year study period (mean absolute change = 20.8 cm +- 18.2 SD). A linear model revealed two contrasting processes of burrow length change: more frequently occupied burrows grew longer, whereas initially longer burrows shortened towards mean burrow length (~50 cm) through collapse. Occupied burrows inside the long-term study site were on average 34% longer and significantly more complex than burrows adjacent to the site, suggesting that researcher disturbance directly or indirectly reshapes burrows. Overall, these results show that storm-petrel burrows are shaped by cycles of renovation and subsequent collapse that may be exacerbated by human disturbance. Given ongoing population declines in storm-petrels and other cavity-nesting seabirds, this study emphasizes the need to balance scientific research with minimizing disturbance, particularly in long-term studies.

Host-microbiota interactions are crucial in adapting to environmental changes, such as fluctuating energy availability. Yet, the extent to which hosts and their microbiota exhibit coordinated plasticity in response to transient energy limitation remains unclear. Here, we examined the capacity for parallel reversible plasticity in hosts and gut microbiota to temporary dietary restriction in hosts with differing energy needs. We subjected bank voles (Clethrionomys glareolus) from a long-term artificial selection experiment for high aerobic metabolism (A-lines) and control lines (C-lines) to dietary restriction (food dilution with insoluble fiber) followed by a recovery period. We measured traits relevant for vole energy consumption and fecal microbiota composition at baseline, after dietary restriction, and after recovery. We hypothesized that A-line voles would experience more severe and prolonged physiological effects of dietary restriction than C-line voles, which might be counteracted by higher microbiota responsiveness. Dietary restriction caused a temporary increase in food consumption, decrease in body mass and resting metabolic rate, and significant shifts in gut microbiota composition, all of which returned to near baseline levels during the recovery period. These effects were largely similar for the selection lines, but contrary to our prediction, the effect of dietary restriction on bacterial diversity was stronger in C-lines. This suggests higher resistance of the A-line microbiota to dietary changes, perhaps resulting from increased host control over microbial communities due to the directional selection for high metabolic capacity. Overall, our findings highlight a substantial capacity for reversible plasticity in host physiology and gut microbiota in fluctuating environments.

Global arthropod decline demands effective biodiversity monitoring strategies. However, most current monitoring approaches do not provide an exhaustive picture of arthropod community structure. In particular, biotic interactions and temporal patterns of biodiversity change are still poorly understood. Here we explore the possibility of addressing these two shortfalls using spiders, one of the most important predators of terrestrial arthropods, as natural samplers for arthropod community DNA. We conducted several experiments comparing the recovered community composition between spider gut contents and traditional monitoring methods. Additionally, we used archived spiders from long-term collections to assess the preservation of prey DNA over time. Spiders proved to be highly efficient natural DNA samplers, with gut content metabarcoding revealing similar community composition and alpha- and beta-diversity like traditional methods. Unique arthropod taxa were detected by spider gut contents and traditional methods respectively, indicating that these two approaches can serve as valuable complements to each other. Well-preserved archived spiders effectively reconstructed historical diets up to 20 years back in time, showing that historical spider collections can act as time capsules of the wider diversity that occurred in the spider's habitat. Spiders as natural samplers overcome critical shortfalls in biodiversity monitoring and contribute to our future understanding of community assembly across space and time.

The toxicity of chemical defences is well studied, but the potential contribution of compound structure to predator deterrence remains largely unexplored. Whether predation acts more strongly on toxicity or unpalatability remains largely untested, partly because few systems allow toxin structure to vary independently of quantity. Heliconius sara larvae provide such a system: those reared on Passiflora auriculata sequester cyclopentenyl cyanogenic glucosides (CGs), while those reared on P. biflora biosynthesise comparable quantities of aliphatic CGs. Using two invertebrate predators, Camponotus floridanus ants and Hierodula membranacea mantids, we tested whether this structural difference affects palatability independent of toxicity. Mantids rejected larvae with cyclopentenyl CGs more often than larvae with aliphatic CGs, despite no detectable difference in total CG content. This pattern was mirrored in extract-based assays with ants, independently of cyanide release: extracts with cyclopentenyl CGs remained deterrent, while extracts with aliphatic CGs did not differ in deterrence from water. Live larvae, by contrast, elicited similar responses from ants regardless of CG structure. These results show that variation in toxin structure can strongly affect palatability, with some compounds conferring greater protection than others. This demonstrates the importance of chemical structural diversity in the evolution of chemical defences.

Recent work on the design of spatial capture-recapture (SCR) studies has largely focused on settings with a single set of movement and detection parameters. In many applied contexts, however, populations are naturally stratified-for example, when different sexes exhibit markedly different movement scales and encounter processes, or when more than one focal species is monitored. Designs that perform well for one group may be inefficient or suboptimal for others, despite the substantial logistical effort typically required to implement an SCR survey. We describe several potential SCR design strategies for surveying multiple groups, and compare these using simulation. Designs considered include regular and optimised grids, clustered and lacework designs, and designs generated using genetic algorithms, including new approaches that optimise information from both groups and a recently published two-stage design. Using proximity detectors such as camera traps, we evaluate design performance in terms of precision, bias, and coverage, and examine how performance changes with the number of available resources. Our results show that multi-group SCR surveys can perform well even when focal groups differ greatly in movement scale and encounter rate, provided that some simple guidelines are followed. Systematic layouts accommodating more than one characteristic spacing e.g. cluster designs, performed consistently well and are a reliable default option. Regular grid designs with a single spacing between detectors are feasible when differences between groups are moderate but performed poorly as differences increased. Establishing when these designs become problematic adds unnecessary risk to the design process. Algorithmic designs did not outperform systematic designs even in highly-skewed scenarios expected to favour them. In practice, when there is more than one focal group we recommend identifying the smallest and largest plausible movement scales and rejecting potentially pathological designs, primarily those that cannot generate recaptures at the smaller scale or cover a substantial area at the larger scale.

Heat-induced sterility is widespread across ectotherm species, yet its underlying mechanisms remain poorly understood. In mated females, thermal stress can disrupt female reproductive physiology while also affecting sperm stored within the reproductive tract, making the mechanisms of fertility loss difficult to distinguish. Here, we investigated the extent and recovery dynamics of female heat-induced fertility loss across seven Drosophila species. We tracked fertility after acute heat exposure, tested whether remating restored fertility, quantified stored sperm in two closely related species, and examined whether reproductive thermal sensitivity covaried with adult heat tolerance. We found that heat exposure reduced egg-hatching success in all seven species, but recovery differed among species. Fertility recovered within several days in most species, whereas D. simulans showed persistent declines and D. melanogaster showed limited recovery. Remating restored or maintained fertility in both species. Stored sperm declined through time in both D. simulans and D. mauritiana, but D. mauritiana retained more sperm at later sampling points, consistent with their contrasting patterns of fertility decline and maintenance. Across species, greater adult heat tolerance was associated with smaller immediate reductions in fertility, but not with cumulative reproductive loss across subsequent days. Together, these findings show that similar initial fertility losses can lead to markedly different recovery trajectories and highlight the importance of comparative, multi-species approaches that integrate female physiology, sperm-storage dynamics, mating behavior, and whole-organism thermal tolerance when assessing reproductive resilience to climate warming.